Ratchet Effect Driven by Coulomb Friction: the Asymmetric Rayleigh Piston 
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The effect of Coulomb friction is studied in the framework of collisional ratchets. It turns out 
that the average drift of these devices can be expressed as the combination of a term related to the 
lack of equipartition between the probe and the surrounding bath, and a term featuring the average 
frictional force. We illustrate this general result in the asymmetric Rayleigh piston, showing how 
Coulomb friction can induce a ratchet effect in a Brownian particle in contact with an equilibrium 
bath. An explicit analytical expression for the average velocity of the piston is obtained in the rare 
collision limit. Numerical simulations support the analytical findings. 
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Introduction - The problem of extracting work from 
unbiased noise is a central issue in the context of en- 
ergy harvesting at small scales [IHIj]. The phenomenon 
of the rectification of non-equilibrium fluctuations is 
called "ratchet effect" , from the seminal works of Smolu- 
chowski Q and Feynman (El, and is studied in the the- 
ory of Brownian motors [6|-|8| . This phenomenon can be 
achieved in the presence of dissipation, i.e. under sta- 
tistical non-equilibrium conditions. It requires that both 
temporal and spatial symmetries are broken. Whereas 
these two fundamental constraints have been pointed 
out the specific mechanisms ruling the action of 
ratchet devices in their different realizations still deserve 
a thorough study within general theories [tl [l(J . 

A class of ratchet models is represented by the "col- 
lisional ratchets" , where fluctuations are induced by the 
interaction of an asymmetric probe with one (or more) 
gas(es) of particles. In these models, realized in several 
variants (lll - [f7l |. dissipation is introduced in two ways: i) 
the probe is in contact with several baths at different tem- 
peratures; ii) interactions between the bath particles and 
the probe are dissipative (granular ratchets). The recti- 
fication of thermal fluctuations has also been observed in 
numerical simulations of an asy mmetric particle diffus- 
ing in a glass forming liquid [18j . In this case, the energy 
flux sustaining the motor effect is induced by the thermal 
unbalance between fast and slow degrees of freedom. 

Here we focus our attention on a different source of dis- 
sipation, introduced by the presence of Coulomb friction 
affecting the dynamics of the probe (or tracer) between 
successive collisions. This kind of dissipation modifies 
the dynamics of the tracer [l9l - |22j ]. and also plays an 
important role in the context of ratchet models [23l428j . 
Here we bring to the fore the structural elements com- 
mon to general collisional ratchets when several forms of 
dissipation are present. We consider a general stochastic 
model ruling the dynamics of the probe's velocity, where 
the interactions with the gas particles are described by a 
Master Equation and the presence of Coulomb friction is 
introduced through a deterministic force. For this model 
we derive an expression for the average drift where two 
terms appear: the first one takes into account the lack 
of equipartition due to non-equilibrium conditions, while 



the second one is directly related to the frictional force. 
The first term represents the heat flux exchanged by the 
probe with the particles bath: it can be induced by a 
non-equilibrium coupling with the bath (e.g. if dissipa- 
tive interactions or reservoirs at different temperatures 
are present), and/or by frictional dissipation. 

To illustrate this result, we study the asymmetric 
Rayleigh piston [2^ in the presence of Coulomb fric- 
tion. In this model a Brownian tracer interacts with 
two gases of particles with different masses but at the 
same temperature. We show that the dissipation through 
friction, coupled with the spatial asymmetry introduced 
by different masses of gas molecules, is sufficient to in- 
duce a net drift on the piston, even if the gases are in 
equilibrium at the same temperature. We stress that 
the motor effect observed in our model originates from a 
mechanism different from those acting in systems where 
the piston is in contact with gases at different tempera- 
tures [T3, H(| , or with ga s particles with different resti- 
tution coefficients [TH, [2(|: in all these systems, indeed, 
non-equilibrium currents are already present even in the 
absence of friction. 

General model- We consider a tracer of mass M and 
velocity V, the motion of which is constrained in one 
dimension, in contact with small particles of one or 
more gases characterized by different parameters: mass 
Wi, temperature Tj, density pi, restitution coefficient 
Ti for the collisions between tracer and gas molecules 
(r.; 6 [0, 1], with Ti = 1 for elastic collisions), where the 
index i denotes different gases. Between successive col- 
lisions with the gas particles the dynamics of the tracer 
is affected by Coulomb friction, slowing its motion. Fur- 
thermore, we assume that the whole system, tracer plus 
gas(es), presents a spatial asymmetry. 

The dynamical evolution of the system is repre- 
sented by a piecewise deterministic process, described by 
stochastic jumps, modeling the interactions of the tracer 
with the particles, and by a deterministic term, taking 
into account the friction. The differential equation de- 
scribing the probability density function of the tracer ve- 
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locity P(V, t) is 



dP(V,t) 
dt 



dV'[W(V\V')P(V',t) - W{V'\V)P(V,t)\ 



A^a(V)P(V,t), 



(1) 



where A is the frictional coefficient, <j(x) is the sign func- 
tion and W(V'|V) are the transition rates for the jump 
from V to V due to collisions, which depend on the mass 
M and the gas parameters (pi, Tj, m, and r{). The spa- 
tial asymmetry may appear in the model through the 
following structure of the transition rates 



W(V'\V) = 



W + (V'\V) 
W-(V'\V) 



if V > V 
if V < V, 



(2) 



with W + (V'\V) ^ W-(V'\V). This structure is realized 
for instance by the probe with triangular shape consid- 
ered in [3, [HI , where collisions on one side always ac- 
celerate the object whereas collisions on the other sides 
slow its motion. The condition (J2J is not sufficient for a 
ratchet effect to be observed. In particular, if A = and 
the transition rates are assumed to satisfy the detailed 
balance (DB) relation with respect to an equilibrium dis- 
tribution Pq(V), then, denoting by (. ..)o the average 
over P Q (V), one has (V) = and (V 2 ) Q = T/M, where 
T is the temperature of the thermal bath of particles (if 
several gases are present, their temperatures have to be 
the same, i.e. Tj = T for every i). 

We stress that model ([1]) is more general than models 
described by Langevin equations (i.e. with white noise), 
widely studied in the literature of Brownian motors 
Indeed, in a Langevin equation non-equilibrium condi- 
tions can only appear through time-varying parameters 
(potentials or temperatures), or through external forces. 
Here, the "kinetic" nature of noise is more physical 
and, due to the different time scales involved, gives the 
possibility to introduce non-equilibrium conditions such 
as constant different temperatures or other dissipation 
channels. 

In this model it is important to put in evidence the 
competition between different timescales: the mean col- 
lision times which, i n the limit of large mass M take the 
form Ti ~ y/rrii/Ti/(piSi), with Si the scattering cross 
section of the tracer with the particles of gas i, and the 
stopping time ta = V* / A due to friction, where V* is 
the average velocity after a collision. These character- 
istic times introduce two regimes in the dynamics. For 
min{Tj} 3> ta, the system is in the rare collision limit, 
where collisions always occur when the tracer is at rest 
because of friction. In this case, the dynamics evolves via 
slip-stick motions and the stationary distribution devel- 
ops a singular contribution in V — 0, as explained in [26| . 
In the opposite limit, when max{r.;} <C ta, the system is 
in the frequent collision limit, and the tracer is never at 
rest. 

Ratchet effect and lack of equipartition- To obtain a 
general expression for the average velocity of the tracer, 



we multiply by V both members of Eq. ([T]) and integrate 
over V. In the stationary state, we get for the momentum 
flow 



= -A(a(V)) + (a(V)}, 



(3) 



where a(V) = f(V- V)W(V'\V)dV is the jump mo- 
ment, which depends on M through the rates W, and 
the symbol (...) denotes an average over the stationary 
distribution P(V). 

For mass of the tracer large enough with respect to the 
largest mass among those of the gas particles, denoted 
hereafter by m, we can do an expansion around V = 
[3l|. Keeping terms up to the second order, we obtain 



a + a'{V) + -a"{V 2 ) 



A(a(V)) ~ 0, 



(4) 



where a' and a" denote the first and second derivatives of 
a with respect to V, respectively, and all coefficients are 
computed in V = 0. In particular, |a/| _1 represents the 
characteristic thermalization time r t h of the tracer with 
the gas, in the absence of friction. This time scale is re- 
lated to the collision time: r t h ~ TiM/rrii [3l|. The coef- 
ficients a, a' and a" are functions of M through the tran- 
sition rates, and have to be expanded in powers of M~ l 
consistently, taking into account that (V 2 ) ~ 0(M~ 1 ). 
Eq. (gj) yields 



(V) 



1 



1 



A 



a(V)) 



= -±[T k -T]+^(a(V)), 
a a 



(5) 



where in the second line we have assumed that thermal 
gradients (if present) are small so that one can define a 
base temperature T, and we have introduced the kinetic 
temperature Tk = M(V 2 ) (assuming (V) 2 <C (V 2 )) and 



a general asymmetry A through the expressions 
a ~ —TA, a"(V 2 )~2AT k . 



(6) 



The above structure for the coefficients a and a" is ver- 
ified in many examples [n], [TH, [HI, O E3 including the 
one discussed below, and follows from Eq. ((SJ [32j (the 
explicit form of A depends on the specific model). 

The interest of Eq. <j5j) is in making clear that there are 
two contributions to the ratchet's drift, corresponding to 
the two channels for heat exchanges of the probe: the 
first one is Dh / = — ^ [Tk — T] , which is proportional to 
the temperature difference T^ — T, and therefore to the 
heat flux exchanged between the ratchet and the thermal 
bath, induced by collisions; the second one is directly re- 
lated to the presence of friction and is proportional to the 
average of the frictional force: Da = Ajp-(a(V)). Notice 
that the first channel can be sustained by the presence 
of reservoirs at different temperatures or dissipative col- 
lisions, but it is also affected by the presence of friction. 
Indeed, if clastic interactions are considered and all the 
baths are at the same temperature, a net flow can still 
be generated by frictional dissipation. 
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FIG. 1. The asymmetric Rayleigh piston with Coulomb fric- 
tion. 



For A = 0, or when the thermalization time r t /j ~ 
1/| a' | is small with respect to the stopping time ta and 
friction can be neglected, only the term D^f remains. 
Then the ratchet effect can be present if and only if the 
transition rates are asymmetric (i.e. A ^ 0) and do not 
satisfy DB, so that the kinetic temperature 71- is dif- 
ferent from that of the external bath T. This is the 
case for many collisional ratchets studied in the litera- 
ture [ll|, US EH EH, E3 j where the explicit expressions 
obtained for the drift in the different cases can be put in 
a form analogous to the first term of Eq. (JSJ) . In the oppo- 
site regime, namely when friction dominates, both chan- 
nels are active and the two contributions may produce 
interesting interplays, with non-monotonic behaviors in 
the drift, as shown below. An expression reproducing 
the non-monotonic drift for collisional ratchets (which is 
a feature already observed, e.g. in (28[) represents a rel- 
evant result of the present study. 

The asymmetric Rayleigh piston- To make more ex- 
plicit our discussion, we now consider the generalized 
Rayleigh piston [29|, [33] , in the presence of Coulomb fric- 
tion. It consists of a piston of mass M, the two faces 
of which are connected with two different gases of elas- 
tic particles of mass m r (at right) and to/ (at left), see 
Fig. [T] The two gases are at equilibrium at the same 
temperature T and have densities pi = p, with i = r.l. 
In such a way the pressures on both sides of the pis- 
ton are equal. The piston velocity is changed by the 
elastic collisions with the (right and left) gas particles 
according to the rule V = V + 1+A // m . (v — V), where 
V and V are the piston velocities before and after the 
collision, respectively. The particles velocities are dis- 
tributed according to the Maxwell-Boltzmann distribu- 
tion Pi(v) = PiV^SF^P ( — ' 2 It _ ) ' where the Boltz- 
mann's constant ks = 1. 

In this model the asymmetric transition rates are [29| 



W+{ V'\V) ( 



(V'-V) 



x pi 



V 2to z 

M + mi , M - mi 



2mi 



-V 



2m, 



V 



, (, > ( M + m r \ 2 , 

w-iy'W) = [-^f) (v-V) 



' M + m r - tr , M — m, 

X Pr I Z V " V 



(7) 



These transition 



2m r 2m r 
rates satisfy DB with re 
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FIG. 2. (color online) DSMC results for (V), with T = 10, 
p — 0.5, m r = 2 and mi = 1, as a function of ta/t^, for 
M = 100 (black dots) and M = 2 (blue squares). The black 
curves show the analytical predictions of the IKM, Eq. (|12|) . 
The dashed lines (red for M = 100 and green for M = 2) 
represent the rhs of Eq. ([5]) with (V 2 ) and (a(V)) computed in 
DSMC. The amplitude of (V) for M = 100 has been magnified 
by a factor 10, for scale reason. Insets: Dhf (red curves) and 
Da (black dashed curves) vs TA/rth- 



spect to the Gaussian distribution Po(V) = 
{2-kT/M)- 1 / 2 exp(-A!V 2 /2T). The explicit expressions 
for the coefficients appearing in Eq. ([5]) are [29l [: 



a = pT[(M + mi)' 1 - (M + m^ 1 ] 
= -pT(mi - m r )/M 2 + C(A/" 3 ), 



a 



-2p 



_M + m t M + m, 
a" = 2p[ mi /(l + mi/M) - m r /(l + m r /M)]/M 



2p(mi - m r )/M + 0{Ar 2 ). 



(8) 
(9) 

(10) 



From Eqs. (jHJ and (|T0|) follows that the explicit formula 
for the asymmetry is A ~ p(mi—m r )/M 2 , which justifies 
the relations ((6j). A similar structure for the coefficients 
a and a" can be traced back in many collisional ratch- 
ets [ll], [TH, [HI, EH, E3] • In this model the time scales are 
ta = V*/A = s/T/M/A, because collisions are elastic, 
and T t h = 1/ 



\a'\ ~ y/ir/{2T)M/[2p(^ml + ^m r 
To study the behavior of the model and to verify the 
relation ([5]) in all regimes, we perform numerical simula- 
tions of the process ([T]) with transition rates (J7J), using a 
Direct Simulation Monte Carlo (DSMC) algorithm [34| . 
We extract the velocity v of a gas particle from pi(v), 
i = r,l with probability 1/2, and then we allow the colli- 
sion with the piston with velocity V to occur with prob- 
ability oc \v — V\. In Fig. [5] (V) is shown (black dots for 
M = 100 and blue squares for M = 2) as a function of 
the ratio T^/rth, which is varied by changing A, with the 
other parameters fixed (see caption) . A net drift is found 
in a wide range of A values: we stress that, at variance 
with kinetic models studied previously, in this system the 
ratchet effect is entirely driven by the Coulomb friction, 
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because the two gases are in equilibrium at the same 
temperature and collisions are elastic. 

The complex non-monotonic behavior of the drift is 
very well described in all the regimes by the r.h.s. of 
Eq. ©, represented in Fig. [2] by red (for M = 100) and 
green (for M = 2) dashed lines. The parameters a, a' 
and a" are given by Eqs. ©, © and (fTU|) and the av- 
erages (V 2 ) and (cr(V)) arc computed in DSMC. The 
behavior of the two terms D^f and Da is reported in 
the insets of Fig. [2] Both terms display plateaux in both 
the opposite limits ta <C and ta 3> t^. The plateau 
in the latter limit is zero for both terms, as equilibration 
with the thermal bath is quickly attained, inducing a zero 
drift. Since also in the opposite limit of rare collisions the 
drift is expected to vanish, as shown below, this produces 
the peaks observed in Fig. [5] Notice that in this model 
Eq. ([5]) also holds for values of M comparable to those 
of the gas particles (see Fig. [5]), if all orders in M are re- 
tained in expressions © , © and ([TU| . This is due to the 
specific forms of the coefficients: in particular, all even 
derivatives of ct(V) greater than the second one vanish 
for this model [29| ■ 

Independent kick model- An analytical explicit formula 
for the average drift can be obtained in the physical sit- 
uation of rare collisions, namely when min{7i} ta- In 
this case, assuming that every collision occurs when the 
piston is at rest, the average velocity can be computed 
in the so-called Independent Kick Model (IKM) (25l.[2r|. 
For our model this yields 



(V) 



dv\v\p r (v) + / dv\v\pi(v 



V(t)dt, (11) 



where V(t) = V - Aa(V a )t, r = \V \/A and V is the 
velocity after a collision: Vo = V + if v > 0, and Vq = V~ 



if v < 0, where V+ = 

l+M jm r 

these expressions one obtains 



and V 



l + M/m, 



{mi + M) 2 (m r + M) 2 



Using 



(12) 



In this formula the net drift explicitly appears when the 
asymmetry in the system is present (i.e. m r ^ mi). For 
small A the formula is not expected to hold because the 
approximation of rare collisions is not valid. Notice also 
that in the limit M — > oo the drift vanishes. In Fig. [2] 
the analytical prediction (|T2|) of the IKM (black lines) 
is shown to be in perfect agreement with the numerical 
results in the rare collision regime. Fig. [5] also shows that 
formula ([5]) is in agreement with the IKM prediction. 

Conclusions- We have presented two interesting re- 
sults: i) formula © for the average drift of a general col- 
lision ratchet in the presence of friction can describe the 
ratchet behavior in all regimes, and explicitly shows the 
two channels of dissipation contributing to the drift; this 
relation has been also tested in a rotor ratchet with dry 
friction recently studied in [28j : ii) Coulomb friction can 
be a source of dissipation sufficient to generate a ratchet 
effect in thermal baths. Our study can be extended to 
other forms of non- linear friction (35l . |36| . 

Our results on the ratchet effect driven by Coulomb 
friction in a thermal bath pave the way to applications 
in the field of nanophysics. At these scales, thermal fluc- 
tuations can be induced by a gas of molecules or a liquid 
environment, and the Coulomb friction is still present, as 
well known from atomic friction experiments 37[ . More- 
over, the developments of new techniques for the design 
and fabrication of nano-devices can provide probes with 
desired shapes and asymmetries. Therefore, all the in- 
gredients are available to realize ratchet devices at small 
scales entirely based on the action of Coulomb friction as 
source of dissipation. 
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